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ABSTRACT 

In this paper we discuss the transport of toroidal magnetic field by a weak meridional flow at the 
base of the convection zone. We utilize the differential rotation and meridional flow model developed 
by Rempel and incorporate feedback of a purely toroidal magnetic field in two ways: directly through 
the Lorentz force (magnetic tension) and indirectly through quenching of the turbulent viscosity, 
which affects the parametrized turbulent angular momentum transport in the model. In the case of 
direct Lorentz force feedback we find that a meridional flow with an amplitude of around 2ms -1 
can transport a magnetic field with a strength of 20 kG to 30 kG. Quenching of turbulent viscosity 
leads to deflection of the meridional flow from the magnetized region and a significant reduction of 
the transport velocity if the magnetic field is above equipartition strength. 
Subject headings: Sun: interior — rotation — magnetic field — dynamo 



1. INTRODUCTION 

Flux-transport dynamos have proven to be suc- 
cessful for modeling the e volution of the large 
scale solar magnetic field (jChoudhuri et alJ 1199.1 
iDikpati fc CharbonneaulU 999(1 . In a flux-transport dy- 
namo the equatorward propagation of the magnetic ac- 
tivity belt (butterfly diagram) is a consequence of the 
equatorward transport of magnetic field at the base of 
the convection zone by the meridional flow. 

However, all studies so far addressed the transport of 
magnetic field by the meridional circulation in a purely 
kinematic regime. The toroidal field strength at the 
base of the solar convection zone inferred from studies of 
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nitude larger than the equipartition field strength esti- 
mated from a meridional flow velocity of a few ms . 
Therefore it is crucial for flux-transport dynamos to ad- 
dress the feedback of the Lorentz force on the meridional 
flow. 

In order to be able to address this question it is nec- 
essary to incorporate a model for the solar differential 
rotation and meridional flow into a dynamo model and 
allow for the feedback of the Lorentz force on differential 
rotation and meridional flow. Differential rotation and 
meridional flow have been addressed in the past through 
mainly t wo approaches: 3D full spherical shell sim- 
ulations llGlatzm aiex-fcGilmanl [1982; Gil man fc Milled 
119861 iMiesch et al.ll2000t iBrun fc Toomrell2002H and ax- 
isymmetric mean field models llKfl chafonoy fc Rucjigerl 
llflM 119951 IRiiHiger et all 119981 IKiiker fc Stixl 12001^ 
While the 3D simulations have trouble reproducing a 
consistent large scale meridional flow pattern (pole- 
ward in the upper half of the c onvection zone) , as 
it is observ ed by helioseismoloev llBraun fc Fan! Il998t 
lHaber et al1l2002t iZhao fc Kosovichevl2004j) . such a flow 
is a common feature in most of the mean field models. 

*The National Center for Atmospheric Research is sponsored by 
the National Science Foundatio n 
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In this paper we present a first step toward building 
a dynamical dynamo model by focusing primarily on 
the transport of a prescribed toroidal magnetic field by 
meridional circulation. We do not attempt in this very 
first model to solve the dynamo equations in order to 
generate the magnetic fields. To this end we will use 
the mean field model for diffe rential rotation and merid- 
ional circulation described bv lRemnell 1(200.11 and extend 
it by incorporating the magnetic tension resulting from 
a purely toroidal magnetic field at the base of the so- 
lar convection zone. In this paper we assume that the 
toroidal magnetic field consists of a homogeneous layer. 
We will briefly discuss in section [21 how an intermittent 
field structure will change the results. In a future paper 
we will present a full axisymmetric mean field dynamo 
simulation. 

2. MODEL 

In this paper we utilize th e differential r otation / 
meridional circulation model of Rempel (2003) and add 
the Lorentz force of a toroidal m agnetic fi eld. F or the 
details of the model we refer to iRemnell 1(200 5). The 
equations we solve are: 
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Pi denotes the pressure perturbation with respect to the 
reference state, po, and si is the entropy perturbation 
s i — Pi/Po — IQi/ffo (made dimensionless with the heat 
capacity c v ), where qq denotes the density of the refer- 
ence state. We use a reference state corresponding to a 



polytropic atmosphere with gravity varying < 



Since 



our model describes the convection zone and overshoot 
region where the deviation from adiabaticity are small 
(|V — V a( j| C 1) we use an adiabatic polytrope for the 
reference state. However small perturbations from adia- 
baticity are considered in the entropy equation Eq. JSJ) 
through the third term on the right hand side. We use 
values for S = V - V ad — -10" 5 below r = 0.725 R Q and 
6 = above. Different profiles of 6 influence primarily 
the differential rotation profile (more specific: the de- 
viation from the Taylor-Proudman state); however, the 
influence on the meridional flow, which is of primary in- 
terest here, is very weak. Qq denotes the core rotation 
rate. n t in Eq. (J5J denotes the turbulent convective heat 
diffusivity and Q the viscous heating. We write the pres- 
sure/buoyancy term in Eq. J5J) assuming small deviations 
from adiabaticity (|V — V a d| < 1). 

F r , Fg, and F$ denote the viscous stress including a 
parameterizati on of turbulent angular momen tum trans- 
port ( A-effect. iKitchatinov & Rudieerl 1(1995^ . The tur- 
bulent angular momentum transport is the driver for 
the differential rotation and the meridional flow in this 
model. The exact form of these terms is discussed in 
detail in lRempell p005|h 

In this model we do not include the induction equation, 
which would be required to address the non-linear evolu- 
tion of the field transported by the meridional flow. As 
a first step towards the full problem we sol ve here for a 
stationary solution similar to lRempel l)2005|h but include 
the magnetic tension force of a toroidal magnetic field. 
This allows determination of the field strength up to 
which a equatorward transport of toroidal magnetic field 
is possible if the feedback of the magnetic tension force 
on the meridional flow is included. For reasons of simplic- 
ity, we omit the magnetic buoyancy term in Eq. (fifth 
term on right-hand side) and focus this discussion only 
on the magnetic tension force. Magnetic buoyancy has 
been a ddre ssed in great detail b y iMoreno-Insertis et al.1 
fl992ft and lRemnel eTaD pOOCh . 

In this study we assume a homogeneous toroidal field 
at the base of the convection zone. Alternatively the 
magnetic field could be highly intermittent leading to the 
following two complications: The Lorentz force is depen- 
dent on the structure of the field in detail and cannot 



be expressed just by the mean field alone; the advec- 
tion velocity of the mean field is not given by the mean 
flow, since field free plasma can flow around magnetic 
elements. In that sense for the case of intermittent field 
there are two extreme scenarios possible: I. The field is 
highly intermittent and couples only weakly through the 
drag force of individual flux elements to the mean flow. 
In this case the equatorward transport of field by the 
meridional flow is mainly dependent on the strength of 
the coupling between mean field and mean flow (similar 
to the behavio r of individual flux tubes as discussed by 
iRemnell ((2003)); the change of the meridional flow pat- 
tern by the feedback of the Lorentz force is a secondary 
effect. 2. The coupling between mean flow and mean field 
is strong. In this case equatorward transport of field by 
the meridional flow is dependent on the change of the 
meridional flow pattern caused by the Lorentz force. 

Even though we assume a homogeneous field in this 
study, our results also apply to intermittent field in the 
case of strong coupling discussed above. Independent 
from the particular way the field couples to the fluid, 
momentum conservation requires that the bulk force act- 
ing on the bulk flow is given by the mean Lorentz force. 
However, the relation between the mean Lorentz force 
and the mean field strength can be more complicated as 
we illustrate in the following very simple example. 

Assume the toroidal field consists of magnetic flux el- 
ements with total filling factor / and field strength Bq. 
The mean tension force (Ft) is proportional to fB 2 , while 
the mean field strength is given by (B) = /-Bo, therefore 
(Ft) ~ (B) 2 j /. This means that for a given tension force 
the assumption of an intermittent field requires a factor 
of f 1 / 2 smaller mean field compared to the homogeneous 
case; however the field strength of individual flux ele- 
ments would be a factor / _1 ^ 2 larger. 

In this paper we evaluate under which conditions an 
equatorward meridional flow at the base of the convec- 
tion zone can transport magnetic field equatorward. To 
this end we compute at a fixed latitude the following ef- 
fective transport velocity for a given magnetic field con- 
figuration: 

^ = 7^. o 

where Vg > means a transport of magnetic flux toward 
the equator. In the case of intermittent field Vg S will be 
an upper estimate for the transport capability. Explicitly 
calculating the coupling between mean flow and mean 
field is beyond the scope of this paper and most likely 
highly dependent on the exact field configuration. 

3. NUMERICAL RESULTS 

In this section we use models for the differential rota- 
tion and mer i dional flow which are very similar to model 
1 in lRempell ll2005[). P arameters that differ from those 
used in iRempell <)2005|) are summarized in Tab. ^ We 
have chosen the parameters such that we have three mod- 
els with meridional flow velocities varying by a factor of 
4, while the turbulent diffusivity is unchanged (model 1, 
2, 3) and three models with a turbulent diffusivity vary- 
ing by a factor of 4, while the meridional flow speed is 
roughly the same (model 1, 4, 5). Fig. 2] shows the flow 
profile of vg at 30° latitude for the models listed in Tab. 
□ 
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TABLE 1 

Parameters of reference models 



model 


Ao 


A [°] 


Kt, v t [m 2 s 1 ] 


1 


0.8 


15 


5 ■ 10 8 


2 


1.49 


7.5 


5 -10 s 


3 


0.43 


22.5 


5 -10 s 


4 


1.36 


7.5 


2.5 ■ 10 8 


5 


0.48 


22.5 


10 9 



Note. — Parameters of differential 
rotation models used. Models 1 to 3 
show different meridional flow velocities 
for the same value of the turbulent diffu- 
sivity, whereas models 1, 4, and 5 show 
roughly the same meridional flow veloc- 
ity for different values of the turbulent 
diffusivity. The parameter Ao deter- 
mines the amplitude of the turbulent an- 
gular momentum flux, the parameter A 
the direction of the angular momentum 
flux with respect to the axis of rotation. 
The quantities Kt and vt denote turbu- 
lent thermal diffusivity an turbulent vis- 
cosity, respectively. 
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Fig. 1. — Meridional flow velocity of the reference models listed 
in Table HI Shown is a radial cut of the velocity at 30° latitude 
of model 1 (solid), model 2 (dotted), model 3 (dashed), model 4 
(dashed-dotted), and model 5 (dashed-triple-dotted). Models 1 - 
3 have the same turbulent viscosity but different flow velocities, 
while models 4 and 5 have similar flow velocities as model 1, but 
different turbulent viscosities. 



Fig. summarizes the results obtained with the ref- 
erence models 1,2, and 3 (varying flow speed, but same 
turbulent viscosity). It turns out that the changes of 
the meridional flow and differential rotation are largely 
independent of the reference model (within a few per- 
cent variation). Therefore we plot in Fig. [21 a) and b) 
the change of meridional flow and differential rotation for 
the reference model 1. In each case we add magnetic field 
with a maximum strength of 1 T, 2 T, 3 T, and 4 T. The 
magnetic field is located at 30° latitude and has a radial 
width of 0.05 Rq and a parabolic profile in B 2 that is 
centered at 0.75 Rq. In latitude we use a Gaussian pro- 
file with a width of about 15°. For the case with 4 T field 
strength the change of meridional flow speed has a maxi- 
mum amplitude of around 7ms -1 and is therefore larger 
than the meridional flow to be expected at the base of 
the convection zone (the observed surface flow speed is 



around 10 — 20 ms -1 , which requires only a return flow 
of 1 — 2 ms -1 at the base of the convection zone due to 
mass conservation). 

A consequence of the meridional flow change induced 
by the magnetic tension force is the formation of a pro- 
grade jet within the magnetized region shown in Fig. 
[2]b), which partially compensates the magnetic tension 
due to the Coriolis force. The amplitude of the jet that 
would be required for a complete compensation is given 
by (see section 0] for a derivation): 

gl = lf "° V (10) 

with the Alfven speed v a = B/y/nop. For a field of 4 T at 
r = 0.75 Rq, 30° latitude, and value of go — 150 kgm -3 
this yields f2^/f2o ~ 0.025, which means that for the 
case shown in Fig [3 the jet compensates about 70% of 
the magnetic tension force (for the other field strengths 
shown, this ratio is about the same). 

The interesting result is that the prograde jet forms 
independent of whether the magnetic field is transported 
equatorward or poleward (angular momentum conserva- 
tion would only provide a prograde jet for a poleward 
movement). This is caused by the fact that the refer- 
ence state (without any influence of magnetic field) is 
characterized by an equilibrium between turbulent an- 
gular momentum transport (including dissipative terms 
and non dissipative terms "A-effect") and the angular 
momentum transport of the meridional flow (the diver- 
gence of the total angular momentum flux has to be zero 
for a stationary solution). The magnetic tension force 
reduces the equatorward flow speed at the base of the 
convection zone, disturbing this balance and therefore 
forcing a change of the differential rotation. Formally 
this result can be understood as follows: Let F° v , 
and be the angular momentum fluxes in the refer- 
ence state due to turbulent non-dissipative angular mo- 
mentum transport (A-effect), turbulent dissipation, and 
meridional flow, respectively. Then the stationary refer- 
ence state is characterized by V • (F A + F° + F m ) = 0. 
The presence of a magnetic field changes the the merid- 
ional flow, leading to a perturbation F m , which changes 
the differential rotation and leads to a perturbation of 
the dissipative angular momentum flux F v . Since the A- 
effect remains unaffected in this case, the new stationary 
equilibrium requires V ■ {F^ + F 1 ) = 0. Since for all 
cases considered F x m is always poleward directed in the 
magnetized region, this requires the formation of a pro- 
grade jet, even though F m +F m can be still equatorward 
d irected. 

iMoreno-Insertis et all ljl992l) : lRemnel et al.ll|2000D dis- 
cussed the formation of prograde jets required for equi- 
librium states of toroidal field at the base of the con- 
vection zone. In their studies a poleward movement of 
the magnetic field was always required, since turbulent 
angular momentum transport and meridional flow were 
not considered in the reference state. The magnetic field 
strength of around 10 T they considered is also signifi- 
cantly larger than the field strength we will focus on in 
this paper. 

While the changes of the meridional flow are indepen- 
dent of the meridional flow speed in the reference model, 
the question whether magnetic field is still transported 
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Fig. 2. — Panel a): Changes of the meridional flow resulting from the presence of a toroidal field with an amplitude of 1 T (solid), 2T 
(dotted), 3T (dashed), and 4T (dashed-dotted) . The solutions shown here use the reference model 1; however, models 2 and 3 show a very 
similar flow pattern (the differences are of the order of a few percent). Panel b) shows the change of differential rotation. The presence 
of a magnetic field leads to the formation of a prograde flow within the magnetized region. The amplitude of this jet is strong enough to 
compensate around 70% of the magnetic tension force through the resulting Coriolis force. Panel c) shows the effective transport velocity 
computed from Eq. J^jJ as function of field strength. Shown are the results for the reference models f (triangles), 2 (asterisks), and 3 
(squares). 
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Fig. 3. — Changes of meridional flow (top) and differential rotation (bottom) using the reference models 4, f, and 5 (left to right). An 
increase of turbulent viscosity (left to right) while keeping the flow speed constant leads to a spread of the side lopes and a reduction of the 
prograde flow within the magnetized region. The linestyle indicates different magnetic field strength as in Fig. [3] The effective transport 
velocity is only marginally affected (see Fig. |5J. 



equatorward depends on the total flow speed. Fig. [21 
c) shows the effective transport velocity according to Eq. 
jjj for the three cases. In the case of a meridional flow of 
around 1.5 ms , magnetic fields up to a about 2.2 T can 
be transported equatorward, while a 6ms" 1 flow could 
transport up to 4.5 T equatorward. Note that the critical 
field strength scales roughly with the square-root of the 
meridional flow velocity. This is different from the case 
of individual flux tubes coupling to the flow through the 
drag force, which results in B ~ v. 
Fig. |21 shows the result for reference models 1, 4, and 



5 that have roughly the same flow velocity, but a vari- 
ation in the turbulent viscosity by a factor of 4. Even 
though the amplitude of the prograde jet changes signif- 
icantly, the transport capability of the meridional flow 
(see Fig. is only marginally affected. In the model 4 
with a value of v t = 2.5 • 10 8 m 2 s _1 the jet compensates 
around 90% of the curvature force (left panels), whereas 
in model 5 with a value of v t = 10 9 m 2 s _1 the jet com- 
pensates only around 50% of the curvature force (right 
panels). This does not affect strongly the transport capa- 
bility of the meridional flow since at the same time with 
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Fig. 4. — Changes of meridional flow (top) and differential rotation (bottom) using the reference models 1 and a decreasing width of the 
magnetic field profile. The width of the field profile centered at r = O.75i?0 decreases from 0.05i?Q to O.O25-R0 (left to right). Similar to 
the increase of turbulent viscosity, a decrease of width leads to a reduction of the prograde flow within the magnetized region. The linestyle 
indicates different magnetic field strength as in Fig. [3] The effective transport velocity is only marginally affected (see Fig. |SJ . 
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Fig. 5. — Effective transport velocity for the cases shown in Fig. 
151 (left 1 and Fig. El fright). Asterisks, triangles, squares correspond 
to the solutions shown in Fig. 151 and Bl on the left, middle, right, 
respectively. 



increasing turbulent viscosity the transport of magnetic 
field through viscous drag increases (a typical value of 
the Reynolds number on scales of the magnetic field is of 
the order of one due to the small meridional flow velocity 
and the large turbulent viscosity). We will discuss this 
in detail in section 

Fig. 0] shows the result for different radial widths of 
the magnetic field. In all cases we use model 1 as the ref- 
erence model. Panels a) and d) show results for a width 
of 0.05 Rq, panels b) and e) for a width of 0.035 R®, and 
panels c) and f) for a width of 0.025 Rq. A comparison 
with Fig. |21 shows clearly that a reduction of the width 
of the magnetic field is equivalent to an increase of the 
viscosity as expected from a simple scaling argument. 

Fig. 13 shows the effective transport velocity for the 



cases shown in Fig. |3 (left panel) and Fig. 01 (right 
panel) . Changes of the turbulent viscosity (by a factor 
of 4) and the width of the magnetic field (by a factor of 
2) considered here influence the critical field strength up 
to which an equatorward transport of flux is possible by 
around 10%. 

4. ANALYTIC ESTIMATE 

Let v r , vg, and f2i be a stationary solution of the differ- 
ential rotation problem without a magnetic field, and v' r , 
v'g, and il'i perturbations around that state caused by the 
presence of the magnetic field. Only considering the 9- 
component of the momentum equation, we can estimate 
the balance between Coriolis force, magnetic tension and 
viscous stress by: 

B 2 v' 

2Q {l' 1 r sin9 cos 9 cot 8- c d z/ t -| = , (II) 

fioQar d 2 

where d denotes a length associated with the radial width 
of the magnetic field and Cd is a coefficient, which takes 
care of the more complicated flow structure. The for- 
mulation given here relates to the more general formu- 
lation of the drag force ~ c w gv 2 /(2d) through the as- 
sumption c w ~ Cd/ Re with the turbulent Reynolds num- 
ber Re — vd/vf Since the large scale flow pattern in 
combination with the large turbulent viscosity is a flow 
with a small turbulent Reynolds number ~ 1, it can be 
expected that the meridional flow behaves like a highly 
viscous fluid. This is justified as long as the scale of the 
flow field is larger than the typical turbulence scale. In 
this limit typical values for Cd should be of the order of 
10 (see text books on fluid dynamics for Stokes law). For 
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small magnetic flux tubes a formulation with c w = const 
rather than c c i = const is more valid. We show later that 
for a reasonable choice of Cd a good agreement between 
this analytic scaling analysis and the numerical result 
can be achieved. 

We emphasize that drawing parallels between laminar, 
viscid laboratory flows and highly turbulent astrophysi- 
cal flows with a low turbulent Reynolds number is specu- 
lative; however it is unavoidable when applying the mean 
field approach, leading to the parameterization of large 
turbulent viscosities. 

The balance between angular momentum transport 
and viscous dissipation yields: 



2Qq — cot I 

r 



^1 „ 



(12) 



We did not introduce here an additional free parameter 
in front of the diffusive loss term, since this prefactor 
~ 1 can be easily absorbed into the definition of the 
length scale d. Combining both equations gives a relation 
between B and v' a : 
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v t dP 



An equatorward transport of the magnetic field requires 
We\ ~ \ v d\i which yields an upper limit for magnetic field 
strength of: 



B ~ /xogo^tan^ (2f2ocos 



The importance of the two terms in the angular brackets 
depends on the value of the turbulent viscosity v t and 
the width of the magnetic field d. For a given width, the 
first term is important for low viscosity, the second one 
for high viscosity. Both are of the same importance if 
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Eq. (fl^jl can be written together with Eq. (fT5|l in the 

form: 
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as suggested 



Eq. i|15|) shows clearly the relation B 
by Fig. 

Inserting Eq. I|13|) into Eq. (|12|l yields for the pertur- 
bation of f2: 



1 
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(17) 



In the case of very low viscosity the perturbation of fl 
is given by 

fil = 1 ( Va 

flo 2 V r2o^ sin ^ 



(18) 



which is exactly the value required to balance the mag- 
netic tension through the Coriolis force. In the case of 
large viscosity the Coriolis force is unimportant and the 
magnetic field is dragged by the meridional flow through 
viscous coupling. 



In order to compare these estimates to the numerical 
results we use the model 1, which is shown in Fig. [3pan- 
els a) and b). As already mentioned earlier the jet ampli- 
tude is sufficient to compensate for 70% of the curvature 
stress, which yields according to Eq. (|17|) v t ~ 0.65^ cr ;t. 
With a value of v t ~ 2.5 ■ 10 8 m 2 s" 1 at 0.75 R Q (roughly 
half the convection zone value) and flo = 2.7 ■ 10~ 6 s" 1 

Eq. (JT5J yields d = cJ^Q.QI^Rq. Using a value of 
Cd = 10 yields d = 0.025 Rq, which is half of assumed 
width of the magnetic field. With Cd = 10, v t ~ 0.65i/ cr ;t, 
go = 150 kgm -3 , r = 0.75 R Q , and v e ~ 2.5 ms -1 Eq. 
(|16|) yields B ~ 2.8 T, which is very close to the numeri- 
cal result. 

Since Eq. (|16|l as function of v t /v cm t has a local min- 
imum for Vt = t'crit, the dependence of B on v t and d 
(through i/ cr i t , Eq. flUO i s expected to be rather weak, 
whereas n^/Oo shows a much stronger dependence as 
found in Figs. and In more detail, Eq. (|17fl docs 
not exactly reflect the scaling of f^/f^o with v t and d 
found in Figs. 01 and ^ which suggests that c<j is depend- 
ing on v t and d itself, however the tendency is indicated 
correctly. 

5. SOLUTIONS WITH QUENCHED VISCOSITY 

So far we have discussed the back-reaction of toroidal 
field on the meridional low through the magnetic cur- 
vature force. Since the equipartition field strength at 
the base of the convection zone is around 1 T (based on 
mixing-length models) the magnetic field also quenches 
the turbulent viscosity to some extent. Since this also 
changes the turbulent angular momentum flux an in- 
fluence on the meridional flow is expected. In order to 
demonstrate this effect use a very weak field of 1 T field 
strength, which has nearly no influence on the meridional 
flow through the tension force, and include a quenching 
of the viscosity given by 

" ;=i "H£>J ■ (i9) 

Note that v t scales in our model the viscous stress and 
the turbulent angular momentum transport (A-effect). 
Fig. [S] a) shows the results for a the values B cq = B max , 
B C q = -Bmax/%/3, and B cq = B max /y/7 that correspond 
to a quenching by a factor of 2, 4, and 8, respectively. 
A comparison with the results presented in Figs. |21to0] 
shows that the back-reaction through quenching of vis- 
cosity has, at least for weak field, a much stronger ef- 
fect than the direct feedback through magnetic tension. 
Whereas in the case of feedback through magnetic ten- 
sion the meridional flow moves around the magnetized 
region on both sides (below and above), in the case of 
quenching the meridional flow closes above the magne- 
tized region. The changes of the differential rotation 
(panel b) are around one order of magnitude lower than 
in Figs. [21 to 01 and no distinct jet forms. 

This results from the fact that the change in the viscos- 
ity changes the parametrized turbulent angular momen- 
tum transport that is the indirect driver for the merid- 
ional flow. In our model the meridional return flow is lo- 
calized in the region where the turbulent viscosity shows 
the strongest radial gradient. The quenching of the vis- 
cosity moves therefore the return flow upward, as clearly 
indicated in Fig. panel c). However, the transport of 
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Fig. 6. — Solution with a quenched viscosity. Panel a) shows the the meridional flow speed, where the dotted line indicates a quenching 
by a factor of 2, the dashed line a quenching by a factor of 4, and the dashed-dotted line a quenching by a factor of 8. The solid line 
indicates the reference solution without any influence of the magnetic field. Panel b) shows the change of differential rotation. Unlike the 
cases shown in Figs. |2*I-IH there is no formation of a significant jet. Panel c) shows the streamlines for the case with a quenching of the 
viscosity by a factor of 8. 



magnetic flux is not switched off completely in this case. 
The magnetic field we assumed is centered at r = 0.75 R@ 
and has a width of 0.05 Rq, which means that the upper 
half is still in a region of considerable flow speed. For ex- 
ample the effective transport velocity is only reduced by 
a factor of two when the turbulent viscosity is quenched 
by a factor of 8). 

6. IMPLICATIONS FOR SOLAR DYNAMO MODELS 

This investigation shows as a very robust result that 
the transport capability of the meridional flow is mainly 
determined by the flow velocity. For a flow velocity of 
the return flow at the base of the convection zone of 
about 2.5ms -1 the maximum field strength that can be 
transported is around 3T (30 kG). This field strength 
is found to be rather insensitive to the turbulent viscos- 
ity or width of the magnetic band. Inspecting Eq. I)l(j[l 
shows that a larger value for this field strength would re- 
quire either a very small or a very large value of v t /f cr it 
or a very large value of a- Since the meridional return 
flow is located roughly where vt shows the strongest ra- 
dial gradient at the base of the convection zone the value 
of i>t in that region will reflect more convection zone val- 
ues rather than very small overshoot values. The value of 
^ cr it is mainly influenced through the effective thickness 
d. Values much larger than the value used in this inves- 
tigation are not feasible, since then the magnetic layer 
would have a larger extent in radius than the meridional 
return flow. For very small values of d an increase in 
the magnetic field strength that can be transported is 
expected. This is not too surprising, since for magnetic 
flux tubes the influence of the drag fo rce is anti propor- 
tional to the diameter. iRempell 1)2003)) showed that flux 
tubes with a diameter of less than 100 km can be trans- 
ported equatorward even if the field strength is around 
10 T, however the magnetic flux associated with these 
flux tubes is orders of magnitude smaller than the flux 
of a typical sunspot. We want to emphasize that for flux 
tubes the use of a ~ Re (c w = const) gives a different 
scaling of B ~ v/d 1 / 2 compared to the B ~ (y t v) l l' 1 jd 
scaling derived from Eq. I)16|) in the limit of small values 
for d. 

Feedback through quenching of turbulent viscosity 



(also affecting the A-effect, which is ~ v t in our model) 
leads to a significant modification of the meridional flow 
if the field strength exceeds equipartition. Since typical 
mixing-length estimates for B e q at the base of the con- 
vection zone are around 1 T, this happens in the same 
field strength range where the direct feedback through 
magnetic tension becomes important, too. Therefore if 
we consider both effects together the effective transport 
velocities are modified, but the results do not change 
dramatically. 

To summarize, the result that the consideration of 
the magnetic curvature force limits the magnetic field 
strength that can be transported towards the equator to 
about 3T (30 kG) seems to be very robust within the 
framework of the mean field model used in this investi- 
gation. 

Current kinematic solar flux-transport dynamos mod- 
els rely on the equatorward transport of toroidal field 
at the base of the convection through the merid- 
ional flow (see compar ison of solutions with and with - 
out meridional flow in [Dikpati fc Charbonneaul l)1999)) : 
Dikp ati"fc Gilmanl 1)200 ll) L The meridional flow veloc- 
ity at the base of the convection zone assumed in most 
of these models is around 1 — 2 ms" 1 , which is close 
to the values considered in this paper. That value is 
also consistent with the observed surface flow of around 
10 — 20ms _1 and a dynamo period around 22 years. The 
main consequence of the work presented here for flux- 
transport dynamos can be summarized as follows: Any 
field exceeding a few T (10 kG) cannot be transported 
toward the equator through a meridional flow with an 
amplitude of a few ms _1 . If the solar dynamo produces 
stronger field (e.g. 100 kG as inferred from studies of 
rising mag netic fl ux tubes llChoudhuri fe Gilmanl Il987t 
Fan et al J 11993 ISchussler et all 11994 ICalieari et all 
1995111998m . the field must get amplified locally through 
induction effects. The s tronger field could be in a n 
equilibrium as discus s ed bvlMor ejio^Jnsertis _et_al] 1)1992)) ; 
IRempel et all 1)20001) : IRempel fc Dikpatil 1)200*37) . which 
would prevent further poleward movement, but an equa- 
torward transport is not possible. 

First attempts to include the feedback on meridional 
flow and differential rotation in a 'dynamic' dynamo 
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model j|Rempel et a l. 2005) show that flux-transport dy- 
namos work with toroidal field strengths up to around 
30 kG, however additional constraints apply through the 
observed limits on the amplitude of tor sional oscillations 
l|Howe et al.ll2004t iRempel et all 12001 . if the feedback 
on differential rotation is included, too. 

The problem addressed in this paper is not limited to 
flux-transport dynamos. Any dynamo will face the prob- 
lem that the toroidal field will start moving toward the 
pole due to the magnetic tension force if the field strength 
is large enough. If the propagation of the magnetic ac- 
tivity is not an advection effect, but a classic dynamo 
wave, the wave will have to compete with the poleward 
movement induced by the magnetic tension force. Since 
the changes of the meridional flow shown in Fig. [5] arc 
only weakly dependent on the meridional flow speed of 
the reference state, they also give an estimate how large 
the poleward movement would be if no meridional flow 
is present. 

We have shown that the formation of a prograde jet is 
an unavoidable consequence when including the magnetic 



tensio n force. Recently iChristensen-D alsgaarc Tet alJ 
( 2004) tried to detect jets associated with toroidal mag- 
netic field in the tachocline. The detection limit they 
found is of the order of 2 — 4 nHz, which is around 
1% of the rotation rate. The jets we see in this study 
have an amplitude of around 1.5% for the strongest 
field we considered (4T). Magnetic field of 2T or less 
leads to the formation of jets with less than 0.5% ampli- 
tude. These jets are therefore at or below the detection 
limit of current helioseismic techniques. If the magnetic 
field has a more complicated intermittent and also non- 
axisymmetric structure it is likely that the amplitude of 
the prograde jet is lower than predicted by our axisym- 
metric model. 



Stimulating discussions and helpful comments on a 
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a very helpful review. 
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